The cyclotron spectrum of anisotropic ultrarelativistic electrons: 
interpretation of X-ray pulsar spectra. 
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The spectrum of cyclotron radiation produced by electrons with a strongly anisotropic velocity dis- 
f>, • tribution is calculated taking into account higher harmonics. The motion of the electrons is assumed 

^-rfH ' to be ultrarelativistic along the magnetic field and nonrelativistic across the field. One characteristic 

^^ feature of the resulting spectrum is that harmonics of various orders are not equally spaced. The phys- 

ical properties and observed spectra of four X-ray pulsars displaying higher cyclotron harmonics are 
analyzed. The cyclotron features in the spectra of all four pulsars can be interpreted only as emission 
lines. Moreover, the observed harmonics are not equidistant, and display certain other properties char- 
(-H ' acteristic of emission by strongly anisotropic ultrarelativistic electrons. In addition, there are indirect 

Oj. theoretical arguments that the electrons giving rise to cyclotron features in the spectra of X-ray pulsars 

I ■ are ultrarelativistic and characterized by strongly anisotropic distributions. As a result, estimates of 

O ' the magnetic fields of X-ray pulsars (which are usually derived from the energies of cyclotron lines) and 

.4.^ , certain other physical parameters require substantial revision. 

1 INTRODUCTION 

^ I A cyclotron line in the spectrum of an X-ray pulsar was first identified for the object Her X-1 in 1977 

il^ ■ [T] . Reports of observations of higher cyclotron harmonics in the spectra of several X-ray pulsars in binary 

^ ; systems (X0115-f63 [5]; Vela X-1 0],|1]; 4U1907-H09 [5], [5]; A0535-h26 [7]j!] were published later. Two 

cyclotron harmonics were also observed in the spectrum of IE 2259+586 [5]; however, the nature of this 
object remains unknown, and will not be discussed here. At least three harmonics have been detected in 
(^ • the spectrum of XO 115-1-63 [9l[T0]. The presence of higher cyclotron harmonics seems to be established 

00 ! most firmly for this source. Unfortunately, the situation for the other objects is less clear, since only weak 

^D ' second harmonics have been detected in their spectra. Nevertheless, even these first observations are very 

interesting, since they can be used to address many theoretical problems connected with the radiation of 
X-ray pulsars. The aim of the present paper is to consider the physical conditions for the formation of 
. , higher cyclotron harmonics in pulsar spectra. We find that a number of problems related to the structure 

j^ ■ of the emitting region of these sources that have remained poorly understood can be easily solved. Here is 

a list of such problems, which will be discussed in more detail below. 

First, it remains unknown if the cyclotron features in the spectra of many X-ray pulsars are absorption 
or emission lines. We will argue below that the four sources indicated above display emission lines (of course, 
if higher-order harmonics are indeed present in their spectra) . 

Second, the velocity distribution of the electrons emitting the cyclotron lines is not entirely clear. We 
can reasonably assume that the velocity of their motion across the magnetic field is fairly low, i.e., weakly 
relativistic (or even nonrelativistic). Otherwise, the cyclotron line would have the characteristic synchrotron 
shape, with a large number of harmonics forming a quasi-continuum. In fact, we can usually see only one 
(the fundamental) harmonic. A few exceptions are listed above, but these sources also possess moderate 
transverse velocities. The electron velocities along the magnetic field are unclear. It is usually assumed 
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that these velocity are small, considerably less than c. Nevertheless, it is possible that the electron-velocity 
distributions are strongly anisotropic, so that the electron velocities are weakly relativistic across the mag- 
netic field and ultrarelativistic along the field. As was shown in the analysis of the spectrum of the X-ray 
pulsar Her X-1 in [llj . there are strong arguments suggesting that the electrons forming the cyclotron line 
are ultrarelativistic, and that their velocities are characterized by very strong anisotropjQ. Unfortunately, if 
only one cyclotron harmonic is observed (as for Her X-1), full verification of this possibility purely through 
spectral analysis is not possible, so that some additional considerations (the dependence of the spectrum on 
the phase of the pulsar, and so on) must be used. The situation is completely different if several harmonics 
are observed. As will be shown below, cyclotron radiation by ultrarelativistic, anisotropic electrons displays 
very specific features. In particular, the harmonics of various orders are not equally spaced. Therefore, by 
analyzing the spectrum, we can determine unambiguously if the motion of the emitting electrons along the 
field is ultrarelativistic or not. 

We will use the following model for the emitting region (hot spot) of an X-ray pulsar, put forward in 
[TT] . When the accretion flow approaches the surface of the neutron star, a shock wave is formed at some 
distance from this surface, and a turbulent region with temperature Tg forms below the shock. The heated 
surface of the neutron star, characterized by the temperature Ts, is located below the turbulent region. The 
ultrarelativistic electrons emitting the cyclotron line are produced at the shock front; their optical depth is 
assumed to be negligible. 

2 CALCULATION OF THE RADIATION SPECTRUM IN A 
COMOVING COORDINATE FRAME 

Let us calculate the magnetobremsstrahlung radiation from electrons located in a constant magnetic field and 
possessing a strongly anisotropic distribution. We assume that the electron velocities along the magnetic- 
field direction (denoted V) are all equal, and that the motion of the electrons along the magnetic field is 
ultrarelativistic. (The situation for the transverse distribution will be specified below.) Let us use a reference 
frame moving along the magnetic field with a constant velocity V . Although this coordinate system is not 
fully associated with the electrons, we will call it the comoving frame, while the initial fixed reference frame 
will be called the laboratory frame. The electrons have only a transverse velocity component v in the 
comoving frame. It is assumed in our model that v <^ c, and the velocity distribution is of the forni^: 

In other words, this is a nonrelativistic, two-dimensional, Maxwellian distribution with temperature T <^ 
nieC^- Let us first find the radiation field in the comoving frame. Further, we will be able to derive the 
distribution of radiation in the laboratory frame by performing a Lorentz transformation. 

According to [12j , the emission of each harmonic in the nonrelativistic case occurs at the fixed frequency 

eH 

LOn — no;// , wheretj// = is the Larmor frequency (2) 

and the intensity of the radiation at the nth harmonic emitted by a particle moving with velocity w <C c 
across the magnetic field 7?, is given by the formula: 
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where 9' is the angle between the magnetic-field vector and the direction toward the observer, J and J 
are the corresponding Bessel function and its derivative, and a prime denotes quantities measured in the 



^Here and below, a strongly anisotropic distribution will be defined as a velocity distribution for electrons whose motion 
along the magnetic field is ultrarelativistic and across the field is nearly nonrelativistic. This subject will be discussed in more 
detail below. 

•^We emphasize that, here and in 1134 l l. v is the total velocity in the comoving frame, not the velocity component perpendicular 
to the magnetic field in the laboratory frame. 



comoving frame. Here and below, the intensity is defined as the energy emitted by the system into unit 
soHd angle per unit time. 

Bessel functions can be expanded into standard Taylor series about zero: 



Thus, if the quantity 
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is not very big, we can take to account only the first term of the expansion. After some transformations we 

obtain: 
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Moreover, in this approximation I 1 — -^ ) — 1, and from ^ we get: 
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where the Thomson scattering cross section ot — —^ ( ^ ) h^s been introduced. To obtain the total 
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intensity of the radiation emitted by N particles distributed in accordance with ([T]), we must integrate ([7]) 
over the distribution: 

oo 
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After dividing this integral by N, we obtain the average intensity of the radiation emitted by one particle: 
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We can see that lim J,'j = oo, when 9' ^ —. Obviously, this is related to the limited applicability of the 

n — *oo 2 

asymptotic expansions of the Bessel function used in ©. In fact, expression ^ represents a zeroth approx- 
imation to the intensity. Taking into account the following terms of the expansion, the condition for the 
applicability of can be formulated as: 

nin + 2) T 

-^-—1 2«0-3 (10) 
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Formula Q yields substantially overestimated values for large transverse temperatures, and is not applicable 
at all when 

!^^«1 (11) 
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As is noted above, the radiation intensities in the laboratory and comoving frames are related by an ordinary 
Lorentz transformation. The comoving frame moves with respect to the laboratory frame with the velocity 
V . Let us introduce the notatioro: 

^ = ^=T? ^=T^h^9 ''=l + Los9' ^'') 

i 7- c c 



■^As in the previous formulas, primed quantities refer to the comoving frame, and unprimed quantities refer to the laboratory 
frame or are frame independent. 



Using the results of [T^, we can easily derive the relations: 
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We also emphasize the important relation 

lift = 7' (15) 

Using p^ - [T5|) . we can easily obtain the relations: 

o 

sind' — — sm9 = — sinOdo' = -^do, (16) 
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where do and do' are elements of solid angle in the laboratory and comoving frames. 

Next, according to |12| . the transformation of the intensity from the comoving to the laboratory frame 
takes the form: 

/ = ^/' (17) 

Radiation at some specific nth harmonic in the comoving frame is monochromatic with frequency hloh In 
the laboratory frame, the frequency of the radiation emitted by one particle is an single-valued function of 
the observation angle 9 (and, therefore, of 9'). Let us find this function in explicit form. According to p^ . 

CO — —uj'. Substituting nuin instead of 9' into this and using (fTS]) . we obtain: 

7 

u) = —nujH = —ntUH = (1 H cos9')jnujH (18) 
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If V^ ~ c and the angles 9' are not very close to tt, we havo 



(I + -COS6I') ~ I + COS6I' (19) 



We obtain in the above approximation 
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sin0' = l-^— -1 =2^—-^— (20) 
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We emphasize that (|20l) is valid only when uj G [0; 2n'yu}H\ (since 9' G [0; tt], so that sa\9' E [0; 1]). 



®The condition for sniallness of the angle, when II19I I. is applicable, can be written: 

(i + cose') - (1 + -^cose') 
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If — ~ 1 this condition is not valid only for angles 9' such that tt — 9' < —. We can see that the solid angle corresponding to 

7 

these plane angles is extremely small ( ~ — J . Since the intensities of all harmonics in the comoving frame have no features 

when 6' ~ TT, the relative fraction of radiation corresponding to this solid angle is obviously also very small j '^ j^ J , and can 
be neglected. 



Since the frequency of the radiation is a single valued function of the angle, we can use this frequency, 
instead of the angle, as the variable specifying the direction. Substituting expression ^ for the intensity 
Z^ into (fTT)) . replacing sin 6*' in accordance with ([20]) . and using the relation ([13]), we finally obtain: 



i^arcH 




(21) 

This formula represents the intensity of the radiation emitted by one particle at the nth harmonic, where 
the directional dependence of the intensity is expressed in terms of the frequency emitted in that direction. 

3 CALCULATION OF THE OBSERVED CYCLOTRON SPEC- 
TRUM 

Let us now find the radiation spectrum detected by an observer an infinite distance from the star. It is 
known (see, for example, [I3])that magnetobremsstrahlung radiation by ultrarelativistic particles is highly 
directional, and is nearly completely concentrated in a narrow cone (with an opening angle of about ~ -) 
along the direction of the particles velocity (in the case under consideration, along the magnetic field). 
Therefore, the radiation detected by an observer at an infinite distance at any particular time is produced 
by a very small part of the emitting region, where the magnetic field is directed precisely toward the observer. 
The linear size of this region is of the order of R-, where R is the stars radius. 

Let us assume that the electrons are uniformly distributed over this part of the emitting region with 
surface density p, that their distribution function has the form Q, and that the magnetic field is constant 
and perpendicular to the stellar surface. Let this surface be separated into bands of width Rd(p, which 
are symmetric with respect to the line from the neutron-star center to the observer. Here, if is the angle 
between rays drawn from the stellar center to the observer and a given point on the stellar surface (i.e., 
the latitude of this point). Let us consider one of these bands. Its area is 2ttR^ sin ipdip, and the number of 
emitting electrons in this band is 

P-2t:R'^ sin ipdip (22) 

Since we have assumed above that the magnetic field is everywhere perpendicular to the stellar surface, the 
angle between the normal to this surface and the direction toward an observer is equal to the angle between 
the magnetic field and the direction toward the observer. As follows from geometrical considerations,: 

ip = e and dip = dO (23) 

The frequency of the radiation emitted by the band is given by (jlSp , where 9 can be replaced by ip: 

w = -nt^H = — y r (24) 

7 7(1-- cos (ys) 

The angle ip is not absolutely constant within the band, and changes by an amount dip. As a result, the 
radiation emitted by the band is not monochromatic, and covers some frequency interval doj. Let us calculate 
the width of this interval. Differentiating ([M]) with respect to (/?, we obtain for — ~ 1: 

duj = uluh — sin ipdip (25) 

7 

Therefore, the band contributes to the total spectrum of the emitting region in the frequency interval 
[w — dui;ti!], where i^ and di^ are given by ([M)) and (|25p . Moreover, it can easily be deduced that the 
corresponding interval for the total spectrum is produced only by electrons in this band. As follows from 
IB 

sm ipdip = doj = (26) 

7 uluh 



Substituting (I26p into (I22p , wc obtain for the number of particles dq emitting in the band: 

dq = 27rpR'^^^ (27) 

/x^ nuJH 

Let do be the sohd angle subtended by the observer at the neutron-star surface. Then, the total energy 
received by the observer from the band under consideration during a time dt will be 

dE = Idqdtdo 

where / is the intensity emitted by one particle. Substituting (^7)) for dq, using p^ . and dividing both 
sides of the equality by dcododt, we obtain an expression for the spectral energy density from the star P„, 
i.e., for the amount of energy emitted by the star per unit frequency interval into unit solid angle per unit 
time: 

P„=/^^^^ (28) 

Dividing this expression by hu;, we make the transformation from the spectral energy density P„ to the 
spectral particle- flux density Qn- Let us introduce the dimcnsionlcss frequency 

(29) 



Substituting formula (I2ip for the cyclotron intensity emitted by one particle into ([28]) and using the relation 

ScarH^ _ ^ _ . L 
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(where A is the fine-structure constant), we finally obtain: 

a. = A,i?^!j(^2^-(^) ) (i+(5-i) M:,:^) (30) 

This formula determines the number of photons of the nth harmonic emitted at some time by the star per 
unit frequency interval into unit solid angle per unit time (i.e., the instantaneous observed spectrum of the 
pulsar). If the physical conditions are the same over the entire surface of the spot, this spectrum coincides 
with the time-averaged spectrum of the pulsar, with a coefficient corresponding to the fraction of time the 
hot spot is observed. Since the transformation (^0]) is valid for (D G [0; 2n], formula (|5n|) is valid when 

< (i < 2n (31) 




Physically, this inequality reflects the fact that, according to (|13p . the photon frequency cannot increase 
by more than the factor 27 in the transformation from the comoving to the laboratory frame, and all 
the photons have the frequency a;„ — uloh in the comoving frame. Therefore, there are no photons with 
frequencies exceeding 2^nu}H in the laboratory harmonic spectrum. Precisely this restriction is responsible 
for the sharp cut-off in the spectrum of the first harmonic. 

We have neglected the gravitational redshift when deriving (pO|) and in subsequent relations. When this 
redshift is fairly small (the stellar radius R is far from the gravitational radius r^), it can easily be taken 
into account by substituting the quantity 



n = Wl--f (32) 

in place of the frequency given by (|29p . Nevertheless, we shall primarily use the previously derived formula 

iOD. 

Combining harmonics of several orders (beginning with the first) and using the constraint (|3ip . we obtain 
the resulting cyclotron spectrum. Here and below, only the first four cyclotron harmonics will be used to 
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Figure 1: Spectrum of cyclotron radiation emitted by an X-ray pulsar containing four harmonics, calculated 
using ((30|) . The transverse temperature of the emitting electrons is T = 20 keV. The dimensionless frequency 
is plotted on a linear scale along the horizontal axis, and the number of particles on a logarithmic scale 
along the vertical axis. 

construct the spectra, while lines of higher orders will be neglected. Let us assume that the distance to 
the pulsar is r = 3.5 kpc, the radius of the neutron star is i? = 10 km (these correspond approximately to 



the parameters of X0115-I-63), and that the surface density of the radiating particles is '--^ 10 



16 



The 



value of h^LUH will be taken to be 10 keV, which is also in reasonable agreement with the energy of the first 
harmonic observed in the spectrum of X0115-I-63 at ^ 20 keV. 

The energy distribution of the photons that should be detected by an observer on the Earth for a 
transverse temperature of the emitting electrons T = 20 keV is presented in Fig. [1] We can see that these 
harmonics are very broad, and overlap each other. In addition, they are clearly not equally spaced; in 
particular, the maximum of the second harmonic nearly coincides with the maximum of the first harmonic. 
The first harmonic possesses a second maximum at uj ~ 0, which is somewhat weaker than the main 
maximum, since higher harmonics are also added to the main maximum. 

It is usual to depict the frequency dependence of the total photon energy, rather than of the number 
of photons. Figure [2] presents this dependence for the same transverse temperature T = 20 keV. The 
formula describing this spectrum can be obtained by multiplying (|30p by hu. As a result, the maxima of 
the harmonics will be shifted; in particular, the maximum of the second harmonic will not coincide with the 
maximum of the first harmonic, and will appear distinct. In addition, the maximum at cj ~ disappears. 
However, as before, the spectral features are not equidistant, and the harmonics are very broad and overlap 
each other. 

These effects have a simple qualitative explanation. Let us consider the emission of photons in the 
comoving frame. First, we must answer the question of which photons will reach the distant observer. It 
is obvious that these will be photons movingaway from the star in the laboratory frame; i.e., those for 
which 6 < ^. According to (fT3)) , this corresponds to the angles < 6*' < tt — i ) in the comoving framqj. 
Therefore, the observer will detect nearly all the photons, apart from those emitted precisely toward the 



'Here and below, we assume that the laboratory frame moves ultrarelativistically with respect to the comoving frame. 



Fig. 2 



1000^ 



100- 



10-1 



> 

CM 

E 
> 

<13 



^ 0,1 

c 




0,01 



2 4 

Dimensionless frequency, co/yoa^ 

Figure 2: Energy spectrum of cyclotron radiation by an X-ray pulsar containing four harmonics. The 
transverse temperature of the emitting electrons is T = 20 keV. The dimensionless frequency is plotted on 
a linear scale along the horizontal axis, and the number of particles on a logarithmic scale along the vertical 
axis. 



star in the comoving frame. 

As was noted above, the energy of all photons corresponding to a single harmonic is the same in the 
comoving frame. Nevertheless, the coefficient of transformation of the photon frequency from the comoving 
to the laboratory frame is not constant, and depends strongly on the angle 0' . Using (fTS]) . we can easily 
show that the frequency of photons emitted precisely toward the observer {9' ~ 0) increases by the factor 
27, the frequency of photons emitted perpendicular to the direction toward the observer (9' ~ 7r/2 in the 

comoving frame) increases by the factor 7, and the frequency of photons emitted at the angle 9' — tt 

decreases by the factor 7. This obviously leads to considerable broadening of the lines. In addition, the 
equidistant character of the lines will be disrupted. This can be interpreted as follows. According to (O), 
most photons associated with the first harmonic are emitted at angles 9' ~ 0. Therefore, their frequency 
increases by the factor 27, and the maximum of the first harmonic occurs at 2jujh- On the other hand, 
most photons associated with the second and higher harmonics are emitted perpendicular to the direction 
toward the observer {9' ~ 7r/2), and particles with 9' ~ are virtually absent (Fig. [3|). The frequencies of 
these photons increase only by the factor 7, so that the maxima of higher harmonics are located at 



ui = 7W„ = jriLUH; n > I 



(33) 



In particular, the maxima of the first and second harmonics coincide. Since the harmonics overlap, their 
maxima are shifted, and relation (|33p is not precisely satisfied. The resulting spectrum differs substantially 
from a spectrum with equidistant features. 

We have considered above only radiation by electrons whose velocities along the field are exactly the 
same. At the same time, a real ensemble of electrons will certainly be characterized by some distribution 
of the momentum component along the field. Therefore, the total velocity distribution of the emitting 
electrons can be written 
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Figure 3: Diagram of the directivity of the first (sohd curve), second (dashed curve), and tenth (dotted 
curve) cyclotron harmonics. The scales are arbitrary. The direction of the magnetic field is marked by the 
arrow. 



This distribution will considerably affect the observed spectrum of an X-ray pulsar. In particular, the first 
harmonic may be transformed from a sharp peak (as in Fig. 2) to a broad line. Therefore, it is interesting 
to answer the following question: can the distribution /(pii) in (|34p be very broad for a real source, and can 
it result in the transformation of separate harmonics into a continuous spectrum? 

As is noted above, ultrarelativistic particles can be produced by collisionless shocks formed in the 
accretion flow. They do not leave the shock, and oscillate at the shock front. Let us suppose that such 
oscillations are roughly harmonic. Then, the particle momentum will vary as p — po cos at, where a 
is the oscillation frequency. Since the oscillations are ultrarelativistic, the momenta of the particles are 
proportional to their Lorentz factors 7 (p ~ 7), so that the 7 values vary in the same way: 



7 = 7o cos at 
The distribution of the corresponding ensemble of particles along the field can be represented 

N 
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(35) 



(36) 



i.e., the distribution is characterized by a sharp peak with its maximum at 7 ~ 79. Of course, real distri- 
butions can differ considerably from (|36p : nevertheless, they can be sufficiently narrow that the individual 
cyclotron harmonics are not transformed into a continuous spectrum. 



4 DISCUSSION 



The parameters of the four pulsars for which higher cyclotron harmonics have been detected are presented 
in the table. Let us consider in more detail the high-energy spectra of these sources, since they contain the 
features of the most interest to us. The spectrum in this region possesses a quasi-power-law character with 
an exponential cut-off (if cyclotron features are not taken into consideration), as is very typical for X-ray 
pulsars. Spectra of this form are supposed to be produced by comptonization of relatively cool radiation 
by hot electrons, considered in detail [13]. The resultingspec trum for energies above the cut-off at E^ 
is approximately a Wien spectrum with the characteristic temperature equal to the temperature of the 
hot electrons Tg [15]. It is formed by photons that have undergone a large number of collisions with the 
electrons, and have therefore been heated to the temperature T^. On the other hand, the photons forming 



Source 


X0115+63 


Vela X-1 


4U1907+09 


A0535+26 


Reference 


[To] 


m 


[5] 


[22] 


Temperature 
of the Planck 
taiire(keV) 


17,4 


9,4 


12 


18,7 


Ratio 

of optical depths 

of the first and second 

harmonics r2/ri 


2,7 


1,8-9,8 


9,3 


2,8 



the power-law part of the spectrum have experienced only a few collisions with hot electrons, and have not 
reached thermal equilibrium. As a result, the power-law spectrum at each specific frequency contains fewer 
photons than a Planck spectrum with temperature Tg. 

To form an absorption line in such a spectrum, the electrons must have a temperature below T^. Even 
if we assume that the electron- velocity distribution differs considerably from maxwellian, the situation is 
unlikely to change fundamentally. In any case, the characteristic energy of motion of the electrons across 
the magnetic field should not exceed T^. The table shows that, for all four sources, the temperature Te 
does not exceed 20 keV. On the other hand, the temperature of the electrons forming complex cyclotron 
lines must be considerably greater: for either emission or absorption harmonics of higher orders to appear 
in the spectrum, the velocity of the electrons across the magnetic field must be comparable to the speed 
of light; i.e., their temperature must be comparable to mc^ ci 500 keV. This requirement can be explained 
quantitatively. The ratio of the optical depths of the ultrarelativistic electrons to cyclotron emission (or 
absorption) at two neighboring harmonics is of the order of the ratio of the mean kinetic energy of motion 
of the particles across the magnetic field to their rest energy. We can easily demonstrate using® that, if 
the motion of the electrons across the field is slow {v <^ c), this ratio is 13 






2n+2 



/ kT 
\2mec'^ 



(37) 



Consequently, if the temperature satisfies the condition 



T <T,r^ 20keV 



(38) 



then the intensity of harmonics should decrease rapidly with increasing harmonic number. In particular, 
the optical-depth ratio for the second and first harmonics should not exceed t- However, observed values 
of this ratio (presented in the table) are 2.7 for X0115-t-63, 9.3 for 4U1907-I-09, > 1.8 for Vela X-1, and 
2.8 for A0535+26 (the cyclotron features were considered in all cases to be absorption features). According 
to ([57)) . this corresponds to electron temperatures of 141, 483, 94, and 146 keV, respectively. It is unclear 
how these hot electrons can absorb radiation with a temperature of ~ 20 keV. Consequently, the cyclotron 
features observed in the spectra of these objects must be emission lines. 

As was reported in [9] and [3] the cyclotron lines in the spectra of XO 11 5+63 and Vela X-1 are not equally 
spaced (the lines were considered in both papers to be absorption features). Let us discuss XO 115+63 in 
more detail. The observed lines in its spectrum [10| are not equidistant (if they are considered to be 
emission features). In addition, the lines are very broad and overlapping. All these characteristic features 
of the cyclotron spectrum are usually explained by the complex geometry of the emitting region (a hot spot 
on the surface of the neutron star). The broadening of the line and its transformation to a broad band are 
thought to be due to variations in the intensity of the magnetic field within the emitting regior|j, while the 



^7In fact, when at least one absorption line is observed in the non-planckian part of the spectrum, formula JSTj is not 
completely accurate, so that the deviations of the source spectrum from a Planck spectrum must be taken into account when 
estimating the temperature of the absorbing electrons from a ratio of relative line intensities. However, for the four objects 
under consideration, our conclusion that the cyclotron features in their spectra are emission features seems to remain valid. 

*As is noted above, in the nonrelativistic case, cyclotron radiation (absorption) occurs at a single frequency, which depends 
only on the magnetic-field intensity, in accordance with l l2l . 
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higher harmonics are formed in parts of the hot spot where the magnetic field is weaker than in the regions 
of formation of the main harmonic, so that their frequencies are shifted with respect to the fundamental 
frequency, iVn < noji. 

Unfortunately, this explanation runs into serious problems. It is not difficult to estimate the relative 
broadening of lines due to the above mechanism. First, we have 

Let us assume that a neutron star is characterized by the standard radius (i? = 10 ) and mass [M = IAMq), 
a dipolar magnetic field 

Hr = ^cose; He^--^sme] (40) 

and a temperature at the foot of the accretion column T^ ~ 10 keV. Then, the height of the hot region is 
approximately determined by the usual barometric formula 

Ar = 



rang 



where k is the Boltzmann constant, nin the mass of the hydrogen atom (the dominant component of the 

accretion gas), and g the free-fall acceleration at the neutron-star surface. Substituting the numerical values, 

we obtain Ar ~ 0.7 . The relative variation in the magnetic field of the form ([^0]) within this distance from 

the neutron-star surface is 

AH Ar 4 

Consequently, ~ 2 • 10~^. Further, let us estimate the variations in the magnetic field over the spot. If 

UJ 

the angular scale of the spot is x, then, as obviously follows from ([iU]) . the corresponding relative variation 
in the magnetic field will be 

AH _ 3 2 

IT ^ 8^ 
According to [T7] , the angular size of the spot is 

X ^ sinx = 7,455 • IQ'^H-'^L^ 

where L is the luminosity of the X-ray pulsar. The luminosity of the pulsar XOl 15+63 is i ~ 10"^^ erg/s 
[2] , while the magnetic field (derived from the energy of the cyclotron line using the nonrelativistic formula) 

would be at least 7 • 10" Gs. Hence, x ~ 0, 06, i.e. — = — — « 1, 4 • 10"^. 

Therefore, the possible broadening or shift of the cyclotron line frequency due to the source geome- 

Aw _, 

try is < 2 • 10 . Of course, this is completely inadequate to explain the observed width of the lines 

UJ 

( ~ 0.2), Moreover, to reproduce the non-equal spacing of the lines, it is necessary to suppose a com- 

plex (and, therefore, quite artificial) temperature distribution over the spot. Thus, we cannot explain the 
observed width of the cyclotron lines, their deviation from an equidistant distribution, or their overlapping 
using only geometrical considerations. The most probable and natural explanation of all these properties is 
that the electrons giving rise to the cyclotron lines are ultrarelativistic, and so have a very anisotropic distri- 
bution of the form (IM)) . Radiation produced by such particles possesses precisely these features. Therefore, 
the presence of several broad cyclotron lines separated by unequal distances in an X-ray pulsars spectrum 
represents a strong argument in favor of our mechanism for the formation of the cyclotron features. The 
spectrum of cyclotron radiation calculated using( ((28| . is plotted in Fig. [4] on a loglog scale for a transverse 
temperature of 20 keV. There is an obvious similarity between the observed spectrum of X0115-I-63 in [lO] 
and the calculated spectrum in Fig. [H Although the first harmonic in Fig. 3] represents a sharp peak, 
whereas it is quite broad in the observed spectrum, this could be explained as an effect of the velocity 
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Fig. 4 
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Figure 4: Energy spectrum of cyclotron radiation by an X-ray pulsar containing four harmonics. The 
transverse temperature of the emitting electrons is T = 20 keV. The scales are logarithmic along both axes. 
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distribution ((34|) of the emitting electrons along the magnetic field. This same distribution should also 
result in additional broadening of the cyclotron lines. 

The calculated cyclotron spectrum displays another very interesting feature: as we can see in Fig. [51 
the dips between neighboring harmonics are located at nearly equal distances from each other. Therefore, 
if these decreases in intensity in the observed spectrum were interpreted as absorption lines, they would 
appear be equidistant. This is probably why the cyclotron features in the spectra of X0115+63, Vela X-1, 
4U1907+09, and A0535+26 have been interpreted as a set of equidistant absorption lines. However, as is 
demonstrated above, they are more likely to be emission rather than absorption lines. 

We have already noted that ultrarelativistic, anisotropic electrons can be produced by a coUisionless 
shock in the accretion flow. As is shown in |15[ I16j such shocks can indeed form, and the corresponding 
electrons will have a strongly anisotropic velocity distribution. This is due to the fact that the electrons 
are accelerated primarily along the magnetic field. In addition, the transverse component of their momenta 
rapidly decreases due to radiative cooling, while the component along the field changes relatively little. The 
estimation of the electron temperature performed in [Kl [IS] yielded values T ~ 10 MeV, corresponding to 
7 ^ 20, i.e. the electrons are ultrarelativistic. 

Note also that there is considerable indirect evidence that the particles forming the cyclotron lines in 
the spectra of X-ray pulsars are ultrarelativistic, primarily associated with discrepancies in estimates of the 
magnetic fields of these objects TT. 

The pulsar Her X-1 was considered in detail in p/7j. In general, its rotation speeds up, but it experiences 
deceleration in some intervals [18]. Consequently, its period is close to the equilibrium period; i.e., the 
angular velocity of rotation of the neutron star together with its magnetosphere is close to keplerian at the 
Alfven radius. We can use this fact to estimate the magnetic field of this pulsar, which turns out to be 
relatively small, about H ^ 3 ■ 10^^ Gs. In addition, radio pulsars in binary systems possess anomalously 
weak magnetic fields, H ~ 10* — 10^^ Gs [11]. Their formation as the result of the evolution of a pair 
containing an X-ray pulsar could be naturally explained if these pulsars had comparatively small magnetic 
fields, H < 3 ■ 10^^ Gs. Mihara et al. [20] explained the observed 35-day cycle of Her X-1 as the result of 
periodic eclipses of the emitting region on the neutron star surface by the accretion disk. For this mechanism 
to operate, the disk must be located at a sufficiently short distance from the surface; i.e., the field must be 
H < 10" Gs. 

At the same time, a cyclotron line at energy 3556 keV is observed in the spectrum of Her X-1. If 
the nonrelativistic formula ([2|), is applied, this corresponds to a magnetic field of (3 — 5) • 10^^ Gs. This 
contradiction appears to be due to the fact that the electrons emitting the cyclotron line are actually 
ultrarelativistic. In this case, as is shown above, the energy of the first (fundamental) harmonic increases by 
the factor 27. Therefore, the nonrelativistic formula gives a substantially (by a factor of 27) overestimated 
value for the magnetic field. Thus, adopting the hypothesis that the velocities of the electrons emitting the 
cyclotron line are ultrarelativistic along the field makes it possible to avoid contradictions between estimates 
of the pulsars magnetic field given by various methods; this represents weighty indirect evidence in favor of 
our mechanism for the formation of the cyclotron radiation. 

Another argument supporting this mechanism is the observation of correlated variations in the energies 
of the cyclotron lines of some pulsars and in their luminosities |21) . This correlation can be explain in a 
natural way by our model. The luminosity of the hot spot is comparable to the Eddington luminosity; 
i.e., the radiation pressure appreciably affects the accretion rate. If the luminosity increases, the accretion 
flow will be substantially slowed. As a result, the intensity of the shock decreases, and the mean energy of 
the ultrarelativistic particles produced by the shock decreases. Consequently, the energy of the cyclotron 
lines also decreases |11| . However, we emphasize that it is recent observations of higher harmonics and 
their deviations from an equidistant spectral distribution that represent the first direct evidence that the 
electrons emitting the cyclotron lines are ultrarelativistic, and have very anisotropic distributions of the 
form ([Mjl. 

In this mechanism for the formation of the cyclotron lines, we must know the electron Lorentz factor 
7 to determine the pulsars magnetic field. In turn, this factor depends on the physical conditions in the 
accretion flow, in particular, in the coUisionless shocks. It is possible that different types of X-ray pulsars 
have different shock structures, and, therefore, different characteristic 7 values. 

Accretion in low-mass binary systems occurs via the flow of material through the inner Lagrange point 
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when the donor star fills its Roche lobe, whereas accretion in massive systems occurs via the capture of 
material from the powerful stellar wind of the OB companion, which does not fill its Roche lobe. It may 
not be a chance coincidence that all four sources displaying several cyclotron harmonics are associated with 
massive systems, so that the second type of accretion is realized. In this case, the conditions at the front of 
the collisionless shock may be more favorable for the ultrarelativistic electrons to acquire greater momentum 
transverse to the magnetic field, resulting in the appearance of higher cyclotron harmonics in the pulsar 
spectrum. 

5 CONCLUSION 

The presence of several cyclotron harmonics in the spectra of X-ray pulsars provides important information 
about their physical properties. 

First, the cyclotron features should apparently not be interpreted as absorption lines. The ratio of the 
optical depths of harmonics of different orders can be used to estimate the temperature of the electrons 
participating in the formation of the cyclotron features. For absorption to take place, this temperature must 
be less than the temperature of the Planck tail that is often observed in the spectra of X-ray pulsars. As 
a rule, Tg is not large: for many pulsars, Tg < 20 KeV. As a result, the second harmonic is substantially 
weaker than the first, and is difficult to detect observationally. In any case, the cyclotron features in all four 
sources for which the second harmonic has been observed are very likely emission features. 

Second, if the cyclotron emission lines are not equidistant, this argues strongly that the electrons emitting 
these lines are ultrarelativistic and have a very anisotropic distribution. Emission by such electrons displays 
characteristic properties that are not typical of emission by nonrelativistic particles, and which are actually 
observed in the observed spectra of X-ray pulsars. There is also other (indirect) evidence supporting this 
mechanism for the cyclotron-line formation. Thus, the cyclotron radiation of X-ray pulsars is most likely 
produced by anisotropic ultrarelativistic electrons. As a result, the nonrelativistic formula ([2|) considerably 
overestimates the pulsar magnetic fields. The magnitude of this overestimation depends on the distribution 
function ([M]) of the emitting electrons in the corresponding shocks. We cannot rule out the possibility 
that this function may be different for different types of X-ray pulsars. Subsequent studies of the physical 
processes occurring in the accretion fiows of these objects may help shed light on this problem. 

Thus, investigations of the spectra of X-ray pulsars containing several cyclotron harmonics may enable 
us to refine estimates of their magnetic fields, and thereby to resolve questions concerning the structure and 
evolution of these objects. 
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